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ABSTRACT 


The  SPINNER  computer  program  has  been  updated  to  compute  aerodynamic 
coefficients  for  a  wide  variety  of  spin  stabilized  projectile  shapes. 
Improvements  over  the  original  program  are  substantial  as  ogive  radius, 
meplat  diameter  and  rotating  band  diameter  are  accounted  for  instead  of 
assuming  mean  values.  Test  cases  are  shown  comparing  the  1969  SPINNER, 
the  1973  SPINNER  and  experimental  data.  Input  instructions  and  sample  pro¬ 
gram  outputs  are  given  along  with  the  1973  program  listing. 
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NOMENCLATURE 


projectile  cross-sectional  area,  ft' 

spin  deceleration  coefficient,  M^/qAd(^j-) 

pitching  moment  coefficient,  M  /qAd 

m 

damping  moment  coefficient,  M  /qAd(qd/2V) 

q_ 

magnus  moment  coefficient,  M  /qAd(pd/2V) 

P 

normal  force  coefficient,  F,,/qA 

N 

magnus  force  coefficient,  F^/qA(pd/2V) 
axial  force  coefficient,  F  /qA 

A 

center  of  gravity,  calibers  from  nose 

2 

axial  moment  of  inertia,  slug-ft 

2 

transverse  moment  of  inertia,  slug-ft 
normal  force,  lbs. 


magnus  force,  lbs. 

axial  force,  lbs. 

spin  damping  moment 

pitching  moment  about  CG 

damping  moment  about  CG 

magnus  moment  about  CG 

total  velocity,  ft/sec. 

projectile  diameter,  ft. 

gravity,  32.174  ft/sec^ 

projectile  mass,  slugs 

projectile  spin  rate,  radians/second 

projectile  pitch  rate,  radians/second 
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NOMENCLATURE  (Continued) 


q 

a 

P 

Subscripts 

a  Derivative  with  respect  to  sin  a 

2 

2  Derivative  with  respect  to  sin 

3 

Derivative  with  respect  to  sin 
oi,.  Derivative  with  respect  to  sin^ 


2  2 

dynamic  pressure,  l/2p  V  ,  lb/ft 
total  angle  of  attack,  radians 

3 

air  density,  slugs/ft 


INTRODUCTION 


The  Armament  Department  of  General  Electric  under  contract  to  Picatinny 
Arsenal  has  developed  an  empirical  computerized  model  for  predicting  the 
aerodynamic  coefficients  of  spin  stabilized  projectiles.  The  code  name  of 
the  new  program  is  SPIN-73. 

The  starting  point  for  the  current  study  was  the  computer  program 
70* 

SPINNER  which  was  developed  at  Picatinny  Arsenal  during  the  period  from 

September  1966  to  October  1968.  This  program  was  modified  by  General  Electric 
68  71 

in  1969  and  1970  to  update  the  predictions  of  the  drag  coefficient  and 
also  to  perform  a  closed  form  dispersion  analysis. 

In  general  the  method  used  during  the  development  of  the  original 
Spinner  was  as  follows:. 

Basic  projectile  configurations  were  selected  which  were  considered  by 
68  70 

Whyte  ’  to  have  well  determined  aerodynamic  coefficients.  Empirical  equa¬ 
tions  and  constants  were  developed,  by  a  trial  and  error  process,  by  which 
the  standard  coefficients  could  be  adjusted  for  changes  in  total  length,  nose 
length,  boattail  length  and  center  of  gravity. 

The  following  limitations  were  and  are  present  in  the  original  program. 

1.  Nose  length  1.8  to  4.0  calibers 

2.  Projectile  length  3.6  to  9.0  calibers 

3.  Boattail  length  0.0  to  1.0  calibers 

4.  Meplot  diameter,  0.10  to  0.15  calibers 

5.  Nose  radius,  secant  +100%  to  secant  -30% 

6.  Rotating  band  diameter,  ] .025  calibers 

*  References  alphabetically  listed  starting  on  page  22. 
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However  as  most  projectiles  in  service  and  under  investigation  during 


the  period  from  1966  to  1969  in  general  fell  within  the  above  bounds  the 
limitations  of  the  program  were  not  considered  very  serious. 

Since  1970  several  programs  have  been  initiated  by  the  Army  and  Navy 
which  are  considering  utilizing  projectiles  with  nose  lengths  of  up  to  5.5 
calibers  and  boattail  lengths  of  up  to  2.5  calibers. 

Also  payload  and  fuzing  capabilities  of  several  small  arms  projectiles 
currently  under  development  by  the  Air  Force,  Navy,  and  Army  have  dictated 
blunter  ogives  (large  meplats)  and  near  tangent  ogives. 

Rotating  band  diameter  are  also  of  larger  scale  on  small  arms  than 
corresponding  shapes  of  large  calibers  thereby  complicating  the  prediction 
process . 

Because  of  these  known  limitations  and  future  requirements  the  need  for 
a  revised  SPINNER  was  indicated.  Thus  this  current  study  was  initiated  in 

August  1972. 

60 

Sears  of  Eglin  in  1972  published  a  computerized  curve  fit  technique 
for  predicting  the  drag  of  projectile.  His  results  indicated  improvement 
over  the  original  SPINNER  in  the  area  of  tangent  ogives  and  meplat  bluntness. 

A  similar  method  to  that  used  by  Sears  was  employed  in  updating  SPINNER. 
In  discussing  the  computer  programs  the  1969  version  of  SPINNER  will  be 
referred  to  as  SPIN-69  and  the  1973  version  as  SPIN-73. 
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PROCEDURE 


The  most  difficult  task  in  the  analysis  of  data  is  determining  a 
constant  accurate  model  which  will  adequately  curve  fit  data  under  all  cir¬ 
cumstances,  whereupon  predictions  of  results  under  a  different  sets  of 
initial  conditions  do  not  result  in  completely  useless  answers. 

An  example  of  useless  results  is  shown  in  figure  15  where  the  pre¬ 
dicted  axial  force  is  negative  in  the  SPIN-69  program.  When  terms  of 

higher  order  polynomials  are  employed  to  obtain  good  fits  one  must  be  very 
cautious  when  using  these  polynomials  to  extrapolate  or  even  interpolate 
data.  These  cautions  are  pointed  out  because  SPIN-73  does  employ  higher 
order  polynomials . 

The  equations  used  for  fitting  and  probable)  errors  v.  xll  be  covered  for 
each  coefficient  individually. 

In  general  the  data  utilized  with  very  fei  exceptions  was  obtained 
from  reports  published  by  the  Ballistic  Resea xh  Laboratory  (BRL)  and 
Arnold  Engineering  Development  Center  (AEDC) .  No  wind  tunnel  data  was  used 

at  all  in  the  data  bank.  Wind  tunnel  data  «ras  used  to  determine  trends 

and  comparisons  were  made  with  the  trends  resulting  data  fitting. 

The  references  utilized  to  collect  the  experimental  data  are  listed 
starting  on  page  22.  Unpublished  data  from  BRL,  Picatinny,  AEDC  and  GE 
were  also  used.  The  method  used  to  curve  fit  the  data  is  described  briefly 
in  Appendix  A. 

Equations  of  the  following  general  form  were  used  for  all  coefficients. 
Definitions  of  VL,  VN,  VB,  VCG,  BD,  DM,  OR,  and  BOOM  are  found  in  figure  1. 
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CX  =  ai  +  a2  (CVN)  +  a3  (CVN  ) 


+  a3  (CXCL)  +  (CXCL  ) 


+  a5  (CVN  •  CXCL) 


+  a,  (CVB)  +  a7  (CBD) 
0  / 


+  ag  (CMK)  +  ag  (CMK  ) 


+  a1Q  (CVN  •  CMK)  +  an  (CXCL  •  CMK) 


+  a..  (CRAT)  +  -  etc. 


where : 

CVN  =  VN  -  2.5 
CXCL  =VL-VN-VB-1.5 
CVB  =  VB 
CBD  =  ED  -  1.02 
CMK  =  CMK  -  1.05 
CRAT  =  VN2/OR  -  0.40 

The  combinations,  variations,  and  parameters  which  can  be  included  in 
the  fitting  equation  are  nearly  infinite.  References  such  as  Dickenson, 

Sears, ^  Watt,^  and  Murphy^  were  used  as  guides  for  determining  the  most 
effective  way  of  deriving  an  empirical  equation.  By  the  trial  and  error  pro¬ 
cess  equations  of  the  above  type  were  manipulated  into  a  form  which  adequately 
described  the  experimental  data. 
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EMPIRICAL  EQUATIONS 


This  section  will  describe  individually  for  each  coefficient  the  equations 
contained  in  the  computer  program  SPIN-73  as  of  June  1973. 


Axial  Force  Coefficient 

CXCL  =  VL  -  VN  -  VB  -  1.5 


CBD 

=  DB  -  1.02 

CDM 

=  (DM  -  0.12) 

2 

CRAT  =  VN2/0R  - 

0.40 

if 

0  < 

VN  <  3.0 

set 

VNX  =  VN 

DXN  =0.0 

if 

0  < 

CXCL  <1.5 

set 

CXCLL  =  CXCL 

.  DXCL  =0.0 

if 

0.2 

<  VB  <  0.65 

set 

VBX  =  VB  -  0.2 

DXBT  =0.0 

If  VN,  CXCL,  or  VB  are  greater  than  the  maximum 
set  VNX  =3.0  ,  DXN  =  (VN  -  3.0)  A^ 

CXCLL  =1.5  ,DXCL  =  (CXCL  -  1.5)  0.01 
VBX  =0.45  ,DXBT  =  (VB  -  0.65)  A1Q 

If  VB  is  less  than  the  minimum 

VBX  =0.0  ,  DXBT  =0.0 

Cx  =  ax  i  a2  (VNX  -  2.5,  +  a3  (VNX  -  2.5)2 

+  a,  (VNX  -  2. 5) 3  +  a  (CXCLL)  +  a,.  (^XCLL)2 

H  JO 

+  a,  (VBX)  +  a„  (CRAT)  +  an  (CRAT)2 

/  O  y 

+  an  (CBD)  +  a  (CDM)  ~  (BOOM/1.36,2  0.01 
-  DXBT  -  DXN  +  DXCL 
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The  "IF"  statements  are  required  to  circumvent  the  need  for  higher  order 
polynomials  in  the  equations.  In  this  manner  only  linear  extrapolation  and 
interpolations  are  allowed  on  the  fringes  of  the  program  capabilities.  This 
should  prevent  completely  erroneous  estimates. 


Normal  Force  Coefficient  Derivative,  Pitching  Moment  Coefficient  Derivative 
and  Normal  Force  Center  of  Pressure 


if  0  <  VN  <  3.0  set  VNX  =  VN  ,  DNX  =0.0 

if  0  <  VB  <  1.0  set  VBNP  =  VBA,  VBMP  =  VBB,  VBX  =  VB 

where:  Subsonic  A  =  1.0,  B  =  0.8 
Supersonic  A  =  1.5,  B  =  1.0 
if  VN  or  VB  are  greater  than  the  maximum 
set  VNX  =  3.0  ,  DNX  =  VN  -  3.0 

VBX  =  1.0  ,  VBNP  =  VB°'5  ,  VBMP  =  VB°’5 


Now  set 

CVNN  =  VNX  -  2.47 

CXLL  =  VL  -  VN  -  VB  -  2.15 

CDMM  =  DM  -  0..17 


CBBD  =  BD  -  1.04 
CCRT  =  VN2/0R  -  0.48 
VBT1  =  CVL/4.7 


CNAB  =  +  B2  (CVNN)  +  B3  (CXLL)  +  B4  (CCRT)  +  B5  (CVNN)2  +  B6  (CXLL)2 

CNBT  =  B?  (VBNP)  +  Bg  (VBX  •  CVNN)  +  Bg  (VBX  •  CXLL) 

CNAT  =  CNAB  +  CNBT 
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AMOMSQ  =  CNAB  +  C2  (CVNN)  +  C3  (CVNN)2 

+  C4  (CVNN)3  +  C5  ( CXLL )  +  C6  (CXLL)2 

+  C?  (CXLL)3  +  Cg  (CCRT)  +  Cg  (CCR T)2 

+  C1Q  (CDMM)  +  Cn  (CCRT  •  CYNN)  +  C1?  (DNX)  ] 

AMOMBT  =  VBTT  [C  2  (VBMP)  +  (VBX  •  CVNN) 

+  C14  (VBX  •  CXLL)  +  C15  (VBX  •  CCRT) 

+  C, ,  (VBX  •  CCRT  •  CVNN) ] 
lb 

CPN  =  (AMOMSQ  +  AMOMBT) /CNAT 

CXT  =  CNAT 
N 

a 

C  ,  =  (VCG  -  CPN)  C.T 

M  N 

a  a 

Yaw  Axial  Force  Coefficient 

CXCL  =  VL  -  VN  -  VB  -  1.5 
CRAT  =  VN2/OR  -  0.40 
CVB  =  VB 

CX  =  D1  +  °2  (CXCL)  +  D3  (crat>  +  d4  (CVB)  ~  CNC( 

The  Yaw  Drag  coefficient  may  be  computed  by  adding  C„  and  C.,  . 

X  N  a 
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CYPA 


CYpct 

Cn  =  (VCG  -  CPP(5)) 
pa(5) 

Damping  Moment  Coefficient 

CLL  =  VL  -  5.0 
CCG  =  VCG  -  3.0 
CVB  =  VB 

C  =  -5.093  [Fn  +  F„  (CLL)  +  F,  (CLL2)  +  F.  (CCG) 

m  i  i  j  4 

9  2 

+  Fc  (CCG)  (CLL)  +  F,  (CCG)  (CLL  )  +  F_  (CCG)  (CVB) 
jo  / 

+  Fg  (CVB) ] 

Spin  Deceleration  Coefficient 

C1  =  G1  (VL/5.51) 

P 

Stability  Analysis 

The  methods  used  for  stability  computations  were  extracted  from  references 
44  and  45.  They  are  identical  to  those  contained  in  the  original  SPINNER. 

Gyroscopic  Stability  Factor, 

2  2,  3  2 

s  =  21  p  /ttI  C  d  V 

g  *  yp  \ 

Dynamic  Stability  Factor,  s^ 

2%  -0  +  <k r2/2)  Cn  ) 
s  _ _ a _ m _ 

d  (C»  -CX-(k2'2/2)C»  +  <kl"2/2)  CJ  ) 

a  q  p 


T7 


. yj&to*  4- 


The  dispersion  (DISP)  is  the  radius  in  mils  of  a  circle  which  a  projectile 
will  impact  in  a  vertical  plane  when  disturbed  to  a  first  maximum  yaw  angle  of 
5  degrees  or  less.  The  basis  for  this  calculation  is  derived  in  Reference  71. 

The  time  step  (DELT)  shown  will  provide  20  integrations  per  nutation 
cycle.  This  is  entirely  adequate  for  a  4th  Order  Range  Kutta  integrator. 
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RESULTS  AND  DISCUSSIONS 


The  results  of  several  test  cases  are  presented  in  Figures  2  thru  15. 
Plotted  are  experimental  points,  SFIN-69  and  SPIN-73  results.  Tabulated 
outputs  of  SPIN-73  are  shown  as  tables  2  thru  15. 

The  following  ranges  of  parameters  are  demonstrated  by  the  test  cases. 


Total  length 
Nose  length 
Boattail  length 
Ogive  radius 
Meplat  diameter 
Band  diameter 


3.8  thru  10.0  calibers 
1.6  thru  5.5  calibers 
0.0  thru  1.0  calibers 
tangent  thru  conical 
0.0  thru  0.26  calibers 
1.00  thru  1.05  calibers 


In  general  the  correlations  between  SPIN-73  and  the  experimental  data 
is  very  good  with  noticeable  improvements  over  SPIN-69.  Most  of  the  effort 
during  this  current  study  has  been  directed  at  the  Axial  Force  and  Pitching 
Moment  correlations  as  these  two  coefficients  are  by  far  the  most  accurately 
determined  during  the  experimental  process.  Much  work  still  remains  to  be 
done  on  these  coefficients  in  terms  of  defining  a  more  adequate  empirical 
model. 

The  most  poorlv  determined  coefficients  remain  the  Magnus  and  damping. 
It  is  this  author’s  opinion  that  the  SPIN-73  improvement  in  these  the  cal¬ 
culations  is  negligible.  While  some  new  data  has  been  published  since  1968, 
in  general  data  was  previously  available  on  similar  shapes.  For  example 
the  projectiles  refered  to  as  the  XM380  and  XM549  were  experimentally 
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investigated  long  ago  as  the  T388  and  T387.  This  data  was  available  in  1967 
and  had  been  included  in  the  original  (SPIN-69)  program. 

The  bulk  of  the  data  published  by  AEDC  through  calendar  year  197?  is 
suspect  as  far  as  the  Magnus  and  damping  coefficient  are  concerned  because 
the  effect  on  linear  theory  reductions  of  a  slowly  varying  pd/2V  was  not  taken 
into  account. 

This  author  also  found  in  several  instances  as  did  Sears  that  the 
geometric  description  of  the  projectiles  under  test  were  not  available 
either  in  the  data  reports  or  the  data  files. 

The  probable  errors  to  the  experimental  data  of  the  SPIN-73  empirical 
equations  are  shown  in  Table  1.  The  number  of  data  points  used  to  compute 
the  probable  ^rror  is  shown  in  parenthesis. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  SPIN-73  program  has  been  shown  through  test  cases  to  be  more  accurate 
than  the  SPIN-69  program. 

The  updating  of  SPIN-73  should  be  continued  as  new  data  is  accumulated. 
Records  should  be  kept  of  shortcomings  and  extremely  poor  predictions. 
Data  should  be  more  carefully  reported  with  respect  to  actual  configura¬ 
tion  tested. 
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CURVE  FIT  TECHNIQUE 


The  method  used  to  perform  least  squares  fits  to  the  experimental 
data  using  the  empirical  equations  was  a  GE  Time  Sharing  computer  program 
code  name  LSQMM.  This  program  was  utilized  with  the  GE415  computer  at 
the  Armament  Systems  Department.  The  brief  description  starting  or  the 
next  page  of  this  program  was  extracted  from  the  following  reference. 

Numerical  Analysis  Routines  //807231A 
Information  Service  Department 
General  Electric  Company 
Bethesda ,  Maryland 

Issued  August  1968  (Revised  Feb.  1969) 


I  J-lL.HIHJJWHW!! 


i!'i.i  hi  ,  s  j 


LSQMM 

T  In',  l  out  mo  determines  tile  coefficients  A(J1,J  =  1,2,  .  . .  N’  of  flic  function 

F(,)-  ■Al,A2Ql,2*  AnQt,n  1  =  U . M 

wfncli  determines  tlie  lies!  appmximation  of  the  function  Y(I)  in  either  the  weighted  least 
squares  sense  or  the  miii-max  sense. 

Usage 

The  Culling  sequence  for  this  routine  is: 

CALL  LSQMMfPi'i  Y  ,  A,RW,M,N,NT,NS,AM,IDIMM,IDIMN) 
where. 


*  PHI  is  the  two  'miensional  array,  PHI(M.N),.  of  coordinate  funcuons  which  are  sup¬ 
plied  by  the  use.  The  k^  column  of  PHI  contains  the  k1*1  coordinate  function  evalua- 
ated  at  each  of  the  data  points  (i.e.,  -  PHI(I.k),  I  -  1,2, M). 

*  Y  is  the  one  dimensional  array,  Y(M),  containing  the  dependent  variables. 

*  A  is  the  one  dimensional  array,  A(N), 

The  A  array  emit  uns  the  loeffieients  A (J)  of  the  function  F(I). 

*  RW  is  the  name 'of  an  array  containing  the  residuals  RfL  Y(I)-F(I),  the  weights  W(I), 
and  temporal  y  storage  to  save  file  vertical  weights  while  doing  the  horizontal  itera¬ 
tions.  It  should  contain  at  least  3*M  locations. 

*  M  is  the  number  of  data  points. 

*  N  is  the  number  of  coefficients,  i.e.,.  number  ol  coordinate  functions., 

*  NT  is  the  maximum  number  of  vertical  iterations.  For  least  squares  fit,.  NT  =  1. 

*  For  least  squares  fit  and  when  there  is  no  division  of  the  data  points  in  min-max,, 
N’S(l)  =  M.  Otherwise,,  NS  is  the  array  containing  the  index  values  of  the  ends  of  the 
sections  when  using  nnn-max  fit .; 

*  AM  is  a  two  dimensional  array,  AM(N,Nj,.  used  internally  to  contain  the  matrix  of 
the  system  of  lm'.ir  equations. 

’  rOIMM  is  th-  first  dimension  of  PHI,  i.e.,  PHLIDIMM.N)'. 

*  IDIMN  is  the  first  dimension  of  AM,,  i.e.,  AM(IDIMN.N). 

Discussion 


If  the  user  wishes  to 
as  follows: 


M 


minimize 


n  1 


1  2 

Y  -F(X  ,  , 

nj  - 


w  '  0,  modify  the  PHI  and  Y 


arrays 


PHI 


y/^T'PHKl.L  ^vTj  PHI(1,2) .  .  y/w,  PHI(1,N) 
^w  PHITOT  ......  .  ...... 


Y 


v'w’A  phhm'o 

V  M 

/yTi] 
y/^TY 72) 
^3_yT3T 


r  wM  PHI(M.N) 
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Sample  Problem 

Find  tii('  second  polynomial  F  A-jX2.  A^v  .  A , which  best  fits  the  iolluwinj' data 

m  the  least  squares  sense,  where  M  -  9,.  N  3,  JIT  1. 

X  -4.,  *3., -2..  -1.,.  0.,.  1.,.  2.:..  3...  4. 

Y  -  2.,  -3.,'  -6  .  -7.;..  -6...  -3.,'  2...  9.,  18. 

Solution  is  F  X'  *  2X-6., 


Sample  Solution 


NEW 

NE4*  FIlE  name  -  -  e  X  ample 
READY  ~  ' 

10  C0MM0N  PH1<9.31,X<9)»A<31»RW(27) 

20  COMMON  AMO, 3), YC9). YA(9).NS< O 
30  DATA  M,N»NT,NS/9,3» 1,9/ 

40  INPUT ,(X(l),l=l«M) 

SO  INPUT, CY<n.l«l»M> 

60  CALL  PH! 1 (M»N ) 

70  CALL  LSQMMCPHI.Y»A»RW,M»N,NT,NS, AM, 9.3) 

80  DO  40  I  =  1 > M 
90  YA<n*0.0 
100  DO  30  J* 1 »N 

110  30  YA(I)=YA<I)+A(J)*PHICI,J) 

120  40  CONTINUE 
130  PRINT  100 

140  100  FORMAT  < '*  X  FCX)  Y-FCX)  AtN) 

150  *  "1 

170  DO  SO  1=1 ,N 

180  SO  PRINT  60»XU),YA<l)»RW<I>,ACI) 

190  60  F0RMAT<4E13.4l 

200  K=N  + 1 

210  DO  70  l*K,M 

220  70  PRINT  80. X< I  >, YA< I 7.RWC 1 > 

230  80  F0RMATC3E13. 4) 

240  STOP 
250  End 

260  SUBROUTINE  PHI  I  CM, N) 

270  C0MM0N  PHI(9,3),X<9> 

280  30  10  1=1, M 

290  10  PHI(I,1)=1.0 

300  IFCN-2)  40,15, IS 

310  15  D0  30  1=1. M 

320  D0  20  J=2 , N 

330  20  PHI C 1, J)=X( I)**< J- 1 • 

340  30  CONTINUE 
3S0  40  RETURN 
■e  EN1 

RUN 

EXAMPI  F 


?'  *  4*  ,  *  >  ~  U«  ,  I  »  •  .  3  *  .  4  * 

72* /  ~  3  • ,"6  *,'7- ,-fc. » -  3  * , 2  * , 9 « ,  18, 


X 

-0.  ’OOOF-t  7, 
-n  ■  t  1 

*0.  2000t  ♦  u  t 

-0. 1000E+01 

o. 

o. 1 000E+01 
0*  20001  +01 
0*3  DOOE *  0 1 
0*400  +01 


FfX!  Y-FCX) 

0* 

JO1,  r*  0  • 

■  60001  -01  0. 

■  -■  •  7  OCOt  ♦  *"*  1  0* 

-0.6000E+01  0* 

-0.3000E+01  0* 

0.2000E+01  0. 

0*  90006  +  01  0* 

0.1R00E+  2  0. 


A  (N  > 

0*6  000E  +  01 
0*  2  000E  +  0 1 
0* 1 OOOE+Ol 
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'JLU 


i.  L  . . (I  I  111  HI  J,! 

* 


INPUT-OUTPUT  DESCRIPTION 

PROCRAM  NAME  -  SPIN-73 


CODING  DATE  -  July  1973 

PURPOSE  -  Predict  the  aerodynamic  coefficients  of  spin  stabilized 
projectiles  at  Mach  numbers  for  0.0  to  5.0. 


Inputs  to 

the  program  are 

the  projectile  physical  dimensions, 

INPUTS 

projectile 
the  local 

mass  properties 
air  temperature. 

,  gun  bore  diameter  and  twist,  and 

Card 

No.  1  (See 

Note  D) 

IBM 

CARD  COL 

VARIABLE 

1 

-  7 

VL 

Projectile  length  -  calibers 

8 

-  14 

VN 

Ogive  length  -  calibers 

15 

-  21 

VB 

Boattail  length  -  calibers 

22 

-  28 

VCG 

Center  of  gravity  -  calibers  for  nose 

29 

-  35 

DM 

bdb 

Diameter  Me'Plat  -  calibers 

36 

-  42 

Rotating  band  diameter  -  calibers 

43 

-  49 

orb 

Ogive  radius  -  calibers 

50 

-  56 

BOOM 

Boom  length  -  calibers 

57 

-  80 

NT1TLE 

Descriptor 

Card 

No.  2  (See 

Note  C) 

1 

-  7 

DIA 

Projectile  diameter  -  inches  2 

8 

-  14 

AX 

Axial  inertia  -  inches  lb-in^  2 

15 

-  21 

TR 

Transverse  inertia  -  inches  lb-in 

22 

-  28 

WOT 

Projectile  weight  -  lbs. 

29 

-  35 

TWIST 

Gun  twist  -  cal/turn 

36 

-  42 

TFftP 

dgunb 

Air  temperature  -  °F 

43 

-  49 

Gun  bore  diameter  -  inches 

50 

-  56 

nautoa 

0  Uses  input,  dimensions 

1  Automatic  dimensions 


A 

BD  set  equal  to  1.02  calibers 
DM  set  equal  to  0.12  calibers 
OR  set  equal  to  2  (-VN2)  (secant  ogive) 

B 

if  input  as  zero  (0,0)  these  inputs  are  changed 
BD  set  equal  to  1.00  calibers 
OR  set  equal  to  secant  ogive 
DGUN  set  equal  to  DIA 

C 

If  aero  estimates  are  only  requirement 
this  card  should  be  left  blank. 

D 

Repeat  cards  1  and  2  to  stack  cases 
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OUTPUTS 


Line 

Line 

Line 

Line 

Line 

Line 

Line 

Line 

Line 


Line 


1 

2 

3 

4 

5 

6 

7 

8 
9 


10 


Organization  designation 
Description  of  item  being  estimated 
Title  line  -  projectile  dimensions 


Projectile  dimensions 

Title  line  -  projectile  physical  properties,  gun  properties 
air  temperature  and  density 

Projectile  physical  properties,  gun  properties,  air  temperature 
and  density 

'Aerodynamic  Coefficients' 

Title  line  -  Mach  No.,  etc. 


Cxo 

CX2 

CNA 

CMA 

CPN 

CYP 

CNPA 

CNPA3 

CNPA5 

CPFl 

CPF  5 

CNPA- 5 


Cmq 

Clp 


-  Zero  yaw  axial  force  coefficient  ^  __ 

-  Yaw  axial  force  coefficient  per  sin  a 

-  Normal  force  coefficient  derivative  per  sin  a 

-  Pitching  moment  coefficient  derivative  per  sin  a 

-  Normal  force  center  of  pressure  -  calibers  for  nose 

-  Magnus  force  coefficient  derivative  per  sin  a 

-  Zero  yaw  Magnus  moment  coefficient  derivative  per  sin  a 

-  Cubic  Magnus  moment  coefficient  derivative  per  sin  a^_ 

-  Quintic  Magnus  moment  coefficient  derivative  per  sin  a 

-  Center  of  pressure  of  Magnus  force  at  1°  yaw  or  less 
calibers  from  nose 

-  Center  of  pressure  of  Magnus  force  at  5°  yaw,  calibers 
from  nose 

-  5°  -  Secant  slope  of  Magnus  moment  coefficient 
derivative  (at  5°  yaw)  per  sin  a. 

-  Damping  moment  coefficient 

-  Spin  deceleration  coefficient 


'Stability  Analysis’ 


GYRO 

SEAR 

RFC  IP 

SBAR5 

RLCIP5 

SPIN 

W1 

W2 

LI 

L2 

LI -5 

L2-5 

DF.LT 

DISP 


Gyroscopic  stability  factor 

Dynamic  stability  factor  at  1°  yaw 

Dynamic  reciprocal  factor  at  1°  yaw 

Dynamic  stability  factor  at  5°  yaw 

Dynamic  reciprocal  factor  at  5°  yaw 

Spin  rate,  radians/second 

Nutation  frequency,  radians/second 

Precession  frequency,  rad  ians/second 

Nutation  damping  factor  per  foot  @  1°  yaw 

Precession  damping  factor  per  foot  @  1°  yaw 

Nutation  damping  factor  per  foot  @5°  yaw 

Precession  damping  factor  per  foot  @  5°  yaw 

Integration  time  step,  seconds  (20  per  nutation) 

Dispersion  factor  per  5°  first  max  yaw,  mils 
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